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Abstract 

A precise determination of the W boson mass at the Fermilab Tevatron requires a theoretical calculation 
in which the effects of the initial-state multiple soft-gluon emission and the final-state photonic correction 
are simultaneously included . Here, we present such a calculation and discuss its prediction on the transverse 
mass distribution of the W boson and the transverse momentum distribution of its decay charged lepton, 
which are the most relevant observables for measuring the W boson mass at hadron colliders. 
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As a fundamental parameter of the Standard Model (SM), the mass of the -boson (Mw) 
is of particular importance. Aside from being an important test of the SM itself, a precision 
measurement of Mw, together with an improved measurement of top quark mass (Mj), provides 
severe indirect bounds on the mass of Higgs boson (Mh)- With a precision of 27 MeV for Mw 
and 2.7 GeV for Mt, which are the target values for Run II of the Fermilab Tevatron collider, Mh 
in the SM can be predicted with an uncertainty of about 35% 111. Comparison of these indirect 
constraints on with the results from direct Higgs boson searches, at the LEP2, the Tevatron 
and the CERN Large Hadron Collider (LHC), will be an important test of the SM. In order to 
have a precision measurement of Mw, the theoretical uncertainties, dominantly coming from the 
transverse momentum of the -boson (P^), the uncertainty in parton distribution function (PDF) 
and the electroweak (EW) radiative corrections to the W boson decay, must be controlled to a 
better accuracy fl. 

At the Tevatron, about ninety percent of the production cross section of W boson is in the small 
transverse momentum region, where < 20 GeV. When Pj^ is much smaller than Mw, every 
soft-gluon emission will induce a large logarithmic contribution to the P^^ distribution so that the 
order-by-order perturbative calculation in the theory of Quantum chromodynamics (QCD) can- 
not accurately describe the P^ spectrum and the contribution from multiple soft-gluon emission, 
which contributes to all orders in the expansion of the strong coupling constant a^, needs to be 
summed to all orders. It has been shown that by applying the renormalization group analysis, the 
multiple soft-gluon radiation effects can be resummed to all orders to predict the P^ distribution 

nil 

that agrees with experimental data ||4l Ij]. RESBOS, a Monte Carlo (MC) program [5] resum- 
ming the initial-state soft-gluon radiations of the hadronically produced lepton pairs through EW 
vector boson production and decay at hadron colliders pp/pp iii2)X, has been used 

by the CDF and D0 Collaborations at the Tevatron to compare with their data in order to de- 
termine Mw- However, RESBOS does not include any higher order EW corrections to describe 
the vector boson decay. The EW radiative correction, in particular the final-state QED correction, 
is crucial for precision measurement of W boson mass at the Tevatron, because photon emission 
from the final-state charged lepton can significantly modify the lepton momentum which is used 
in the determination of Mw In the CDF Run lb W mass measurement, the mass shifts due to 
radiative effects were estimated to be —65 ± 20 MeV and —168 ± 10 MeV for the electron and 
muon channels, respectively fl. The full next-to-leading order (NLO) 0{a) EW corrections have 
been calculated and resulted in WGRAD ||2|], a MC program for calculating 0{a) EW ra- 



diative corrections to the process pp i'ei{'y). However, WGRAD does not include the dominant 
correction originated from the initial-state multiple soft-gluon emission. To incorporate both the 
initial-state QCD and and final-state QED corrections into a parton level MC program is urgently 
required to reduce the theoretical uncertainties in interpreting the experimental data at the Teva- 
tron. It was shown in Refs. 1^1^ that at the NLO, the EW radiative correction in pp ii^ii'y) 
is dominated by the final-state QED (FQED) correction. Hence, in this paper we present a con- 
sistent calculation which includes both the initial-state multiple soft-gluon QCD resummation and 
the final-state NLO OED corrections, and develop an upgraded version of the RESBOS program, 
called RESBOS-A llO], to simulate the signal events. 

The fully differential cross section for the production and decay of the W boson that includes 
only the effect of the initial- state multiple QCD soft-gluon emission can be found in Ref. jsl. 
To include also the final-state NLO QED contributions, we sum up the following two sets of 
differential cross sections. One, cf. Eq. ([T]), contains final-state QED virtual correction and part 
of the real photon emission contribution in which photon is either soft or coUinear. Another, cf. 
Eq. Q, includes the hard photon contribution from the real photon emission processes. Denote Q, 
y, Qt and (pw to be the invariant mass, rapidity, transverse momentum and azimuthal angle of the 
di-lepton pair, respectively. For production and decay, we have 



+ F(gr,Q,Xi,X2,C4)|, (1) 

and 

dQ^dydQld<Pwdn^ J,,, 

where and dlla represents the two-body and three-body phase space of the vector boson 
decay products, respectively. In the above equations the parton momentum fractions are defined 
as xi = — ^ and X2 = — where is the center-of-mass energy of the hadrons hi and /12. 
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The renormalization group invariant quantities Wj'^\b) and W^'^\b), which sum to all orders in 
as all the singular terms that behave as Q^"^ x [1 or \n(Q'^/Q'^)] for Qt 0, are 
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and 
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where (Xp is the Born level parton cross section for jk vd'^, and crp''^^ includes the final-state 
NLO QED corrections. The notation ® denotes the convolution jsl 

and the V^^ coefficients are the Cabibbo-Kobayashi-Maskawa mixing matrix elements. In the 
above expression j represents quark flavors and k stands for antiquark flavors. The indices a and 
b are meant to sum over quarks and antiquarks or gluons. Summation on these double indices 
is implied. As compared to the results shown in Refs. flQ], W^^{b) and w')^{b) contain addi- 
tional ^ corrections, which come from the final state QED corrections. The Sudakov exponent 
S{b, Q, Ci, C2) in Eqs. © and © is defined as |3] 
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The explicit forms of the A, B, and Cjn functions and the renormalization constants Ci (i = 1, 2, 3) 
can be found in Appendix D of Ref. o]. In our calculation, we have included A^-^^ , A*^^) , i?*^^) , 5*^^^ , 
C'^^^and C''^\ with canonical choice of Cj's. The Y piece in Eq. ([l]), which is the difference of the 
fixed order perturbative result and their singular part, can be found in Appendix E of Ref. jsl. 

We follow the prescription in Ref. |7], which decomposes the electroweak 0(a) contribution 
to the resonant single W production in a general 4-fermion process into gauge invariant QED-like 



and weak parts, to extract a gauge invariant QED-like form factor from the photon contribution. 
The NLO FQED differential cross sections are calculated by using phase space slicing method lO], 
which introduces two theoretical cutoff parameters, soft cutoff Ss and coUinear cutoff Sc, to iso- 
late the soft and coUinear singularities associated with the real photon emission subprocesses by 
partitioning phase space into soft, coUinear and hard regions such that 

The soft region is thus defined by requiring that the photon energy (E^) in the (jk) parton center- 
of-mass frame to satisfy < 5sV§/2, where is the invariant mass of the (jk) partons. Using 
the dimensional regularization scheme, we can then evaluate, in n-dimensions, the real photon 
emission diagrams under the soft photon approximation, where the photon momentum is set to be 
zero in the numerator, and integrate over the soft region. In the soft and coUinear regions the cross 
section is proportional to the Born cross section. The soft singularities originating from the final- 
state photon radiation cancel against the corresponding singularities from the final-state virtual 
corrections and leave a finite result depending on the soft cutoff parameter 6s- For E^ > 6sV^/2, 
the real photon emission diagrams are calculated in four dimensions using the helicity amplitude 
method. The coUinear singularities associated with photon radiation from the final-state charged 
lepton is regulated by the finite lepton mass. The end result of the calculation consists of two sets 
of weighted events corresponding to the jk ^ ue£~^('j) and jk ^ i^d'^'j contributions which are 
included in Eqs. ^ and Q, separately. Each set depends on the soft cutoff parameter 5s. The sum 
of the two contributions, however, is independent of 5s, as long as the soft cutoff is small enough 
to validate the soft-gluon approximation. In our numerical studies, we take 5s = 0.001 which 
yields a stable numerical result in agreement with Refs. 1^0]. Through our calculation, we adopt 
the CERN LEP line-shape prescription of a resonance state and write the W boson propagator as 

where Tyy is the width of W boson. 

To examine how much the combined contributions from the initial- state QCD resummation and 
the final-state QED corrections can affect the precision measurement of Mw, we perform Monte 
Carlo analyses to study a few experimental observables that are most sensitive to the measurement 
of Mw at the Tevatron (a pp collider with ^/S=\.96 TeV). For the numerical evaluation we chose 
the following set of SM parameters: a = 1/137.0359895, = 1.16637 x lO^^GeV-^, Mw = 
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FIG. 1: Transverse mass distribution of boson 

80.35 GeV, Vw = 2.0887 GeV, Mz = 91.1867 GeV, m, = 0.51099907 MeV. Thus, the square of 
the weak gauge coupling is g'^ = Ay/2M^,G^. Because of the limited space, we focus our attention 
on the positively charged electron lepton (i.e. positron) only, though our analysis procedure also 
applies to the /i lepton. The complete study including both electron and muon leptons will be 
shown in our forthcoming paper, in which we also extend our study to the LHC. 

The W events in these analyses are selected by demanding a single isolated high px charged 
lepton in conjunction with large missing transverse energy. To model the acceptance cuts used by 
the CDF and D0 Collaborations in their W mass analyses, we impose the following transverse 
momentum (jjt) and pseudo-rapidity (rj) cuts on the final-state leptons: 



> 25GeV, |7;(e)| < 1.2, > 25GeV. 



(9) 



Due to the overwhelming QCD backgrounds, the measurement of W boson mass at hadron 
collider is performed in the leptonic decay channels. Since the longitudinal momentum of the 
neutrinos produced in the leptonic W boson decays (W'^ i^e^'^) cannot be measured, there is 
insufficient information to reconstruct the invariant mass of the W boson. Instead, the transverse 
mass distribution of the final state lepton pair, which exhibits a Jacobian edge at Mt ~ M^, is 



used to extract out Mw- Transverse mass (M^) of W is defined as 



= -y2p^p^(l - COS0) , (10) 

where (p is the angle between the charged lepton and the neutrino in the transverse plane. The 
neutrino transverse momentum (p^) is identified with the missing transverse energy (^t) in the 
event. In Fig. [U we show various theory predictions on the M]f distribution. The legend of the 
figure is defined as follows: 

• LO : including only the Bom level initial-state contribution, 

• RES : including the initial-state multiple soft-gluon corrections via QCD resummation, 

• LO QED : including only the Bom level final-state contribution, 

• NLO QED : including the final-state NLO QED corrections. 

For example, the solid curve (labelled as RES-i-NLO QED) in Fig. [Ha) is the prediction from our 
combined calculation given by Eqs. ^ and @. 

As shown in Fig. da), compared to the lowest order cross section (dotted curve), the initial 
state QCD resummation effects (dashed curve) increase the cross section at the peak of the 
distribution by about 20%, and the final state NLO QED corrections (dot-dashed curve) decrease 
it by about —12%, while the combined contributions (solid curve) of the QCD resummation and 
FQED corrections increase it by 7%. In addition to the change in magnitude, the line- shape of 
the M]f distribution is significantly modified by the effects of QCD resummation and FQED 
corrections. To illustrate this point, we plot the ratio of the (RES-i-NLO QED) differential cross 
sections to the LO ones as the solid curve in Fig. db). The dashed curve is for the ratio of 
(LO-i-NLO QED) to LO. As shown in the figure, the QCD resummation effect dominates the 
shape of M]f distribution for 65 GeV < Mw < 95 GeV, while the FQED correction reaches its 
maximal effect around the Jacobian peak (Mjf ~ Mw)- Hence, both corrections must be included 
to accurately predict the distribution of around the Jacobian region to determine Mw We 
note that after including the effect due to the finite resolution of the detector (for identifying an 
isolated electron or muon), the size of the FQED correction is largely reduced \ ^^. 

Although the M^ distribution has been the optimal observable for determining Mw at the 
Tevatron, it requires an accurate measurement of the missing transverse momentum direction 
which is in practice difficult to control. On the other hand, the transverse momentum of the decay 
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FIG. 2: Transverse momentum distributions of 

charged lepton (pf,) is less sensitive to the detector resolution, so that it can be used to measure 
Mw and provide important cross-check on the result derived from the Mjf distribution, for they 
have different systematic uncertainties. Another important feature of this observable is that 
distribution is more sensitive to the transverse momentum of W boson. Hence, the QCD soft- 
gluon resummation effects, the major source of , must be included to reduce the theoretical 
uncertainty of this method. In Fig. 12 a), we show the pf. distributions predicted by various theory 
calculations, and in Fig. E^b), the ratios of the higher order to lowest order cross sections as a 
function of pf.. The lowest order distribution (dotted curve) shows a clear and sharp Jacobian peak 
at pfp ~ Mw/2, and the distribution with the NLO final-state QED correction (dot-dashed curve) 
also exhibits the similar Jacobian peak with the peak magnitude reduced by about 15%. But the 
clear and sharp Jacobian peak of the lowest order and NLO FQED distributions (in which = 0) 
are strongly smeared by the finite transverse momentum of the W boson induced by multiple soft- 
gluon radiation, as clearly demonstrated by the QCD resummation distribution (dashed curve) and 
the combined contributions of the QCD resummation and FQED corrections (solid curve). Sim- 
ilar to the M]f distribution, the QCD resummation effect dominates the whole range, while 
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FIG. 3: Distributions of TZ, Si and 62 



the FQED correction reaches it maximum around the Jacobian peak (half of Mw)- The combined 
contribution of the QCD resummation and FQED corrections reaches the order of 45% near the 
Jacobian peak. Hence, these lead us to conclude that the QCD resummation effects are crucial in 
the measurement of from fitting the Jacobian kinematical edge of the distribution. 

It is also interesting to examine the effect of the final-state NLO QED correction to the the- 
ory predictions with LO (leading order) or RES (resummed) initial-state cross sections, which is 
described by the observable TZ, defined as 
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where 
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The distributions of 7^ as a function of and M)f are shown in the upper part and lower part of 



Fig-El respectively. As expected, the 1Z{Mt) distribution is almost flat all the way from 60 GeV 
to 100 GeV, but the TZ{pt) distribution deviates from one in the Jacobian peak region. This is due 
to the fact that is more sensitive to Pj^ than . 

In order to study the impact of the presented calculation to the determination of the W boson 
mass, the effect due to the finite resolution of the detector should be included, which will be 
presented elsewhere. 
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